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ABSTRACT 

The physical characteristics of the Lyman-a forest cloud systems are com- 
bined with observations on the baryonic mass density of the Universe and con- 
straints from primordial nucleosynthesis to set boundary conditions on the In- 
tergalactic Medium (IGM) at the epoch oi z = 2.5. The Universe is consid- 
ered a closed system and allowed to expand adiabatically from the epoch when 
QSOs first ionized the IGM ( 5 < Zon < 20). The average kinetic energy of 
a gas is calculated in the region where the gas transitions from relativistic to 
non-relativistic behavior. All of the above measurements are then used to de- 
termine the thermal history of the IGM in the redshift range 2.5 < z < Zon- 
The hot IGM is assumed to inverse Compton scatter photons from the Cosmic 
Microwave Background (CMBR) and consequently distort the CMBR as seen 
at the present epoch. The temperature of the IGM at z = 2.5 and the epoch 
Zon are adjusted, within the constraints defined above, to give the best overall 
agreement with published data on the temperature of the IGM. We find that the 
model of the IGM proposed here does not grossly distort the CMBR, and in fact 
agrees quite closely with the preliminary results from the Cosmic Background 
Explorer (COBE) satellite. However, our model of the IGM cannot explain the 
observed cosmic x ray background. 

This paper was originally written in 1990. It was never submitted for publi- 
cation. 



1. Introduction 



In the late 1940s a "Big Bang" theory of nucleosynthesis was postulated by George Gamow 
and his students, Ralph Alpher and Robert Herman. In 1948, Alpher and Herman used 
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this theory to predict a blackbody background with a present temperature of about 5 K. 
Similar calculations were done with similar results by Ya. B. Zeldovich, Fred Hoyle and 
R.J. Taylor in 1964. Independent of the above calculations, P.J.E. Peebles argued in 1965 
that there ought to be a blackbody background of radio noise left over from the early 
Universe, red shifted to a present temperature of about 10 K. Peebles based his argument 
on the fact that if there had not been an intense radiation field during the initial few 
minutes of the expanding Universe, that nucleosynthesis would have proceeded at such 
a rate that a substantial amount of the primordial hydrogen would have been fused into 
heavier elements. This would be in direct conflict with the fact that roughly 75 percent of 
the observed baryonic matter of the Universe is hydrogen. The intense radiation field was 
needed in order to prevent the hydrogen nuclei from fusing. 

This radiation field would have survived, albeit at lower and lower temperatures as the 
Universe continued to expand. Therefore this radiation would be present today and if there 
had been no interaction of the radiation field with other constituents of the Universe, the 
radiation field would be isotropic, homogeneous and a perfect blackbody. In 1964, Arno 
Penzias and Robert Wilson discovered this relic radiation field using a 20 foot horn reflector 
antenna located on Crawford Hill at Holmdale, New Jersey. 

Since 1965 a variety of experiments have been performed in order to confirm that Penzias 
and Wilson actually observed the relic radiation of the Big Bang. It is hypothesized that 
this relic radiation, with a temperature of T^o ~ 3 ET, would be dominate in the wavelength 
range 1 cm > A > 0.3 mm. The wavelength at which the maximum intensity would occur 
is found from the Wien displacement law 

Amax7; = 0.29 . (1) 

For = 3 ii', Xmax ~ 1 mm. Unfortunately this entire range of wavelengths is not 
observable from the earth's surface. Molecular absorption bands in the earth's atmosphere 
prevent observations at wavelengths less than 1 cm. Rocket borne instruments are needed 
to observe the Wien Tail of the relic radiation. 

A number of groups have measured the intensity of the relic radiation and have found it to 
be remarkably isotropic (discounting the dipole moment observed because of our galaxy's 
peculiar motion in the direction of the Virgo cluster of galaxies), and well fit by a black- 
body except in the Wien Tail of the spectrum, A < 1 mm. This data is shown in Table [T] 
below. As seen in Tabled! measurements of T^, with the exception of the COBE data, indi- 
cate a considerable short-wavelength excess in the cosmic microwave background radiation, 
CMBR. Many attempts were made to explain blackbody distortions of this type, including: 

i. Very early superstars or hypercactive protogalaxies released extraordinary amounts 
of radiant energy that was then reradiated in the infrared by a prodigious density of 
cosmic dust. 
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Wavelength (cm) 


TcMBR {K) 


Reference 


50.0 


2.45 ± 0.70 


Sironi et al. (1987) 


21.2 


2.28 ±0.39 


Levin et al. (1988) 


12.0 


2.79 ±0.15 


Sironi & Bonelli (1986) 


8.1 


2.58 ±0.13 


DeAmici et al. (1988) 


6.3 


2.70 ±0.07 


Mandolesi et al. (1986) 


3.0 


2.61 ±0.07 


Kogut et al. (1988) 


0.909 


2.81 ±0.12 


Smoot et al. (1985) 


0.333 


2.60 ±0.10 


Johnson & Wilkinson (1987) 


1.2 


2.783 ±0.025 


Crane et al. (1989) 


0.264 


2 796+°-°^^ 


Crane et al. (1989) 


0.132 


9 7c;+0.24 
— u.zy 


Crane et al. (1986) 


0.264 


2.70 ± 0.04 


Meyer & Jura (1985) 


0.22 


2.60 ± 0.09 


Bersanelh ei al. (1989) 


0.132 


2.76 ±0.20 


Bersanelh ei aZ. (1989) 


0.116 


2.799 ±0.018 


Matsumoto et al. (1988) 


0.071 


2.963 ±0.017 


Matsumoto et al. (1988) 


0.048 


3.150 ±0.026 


Matsumoto et al. (1988) 


0.351 


2.80 ±0.16 


Peterson et al. (1985) 


0.198 




Peterson et a/. (1985) 


0.148 


2.92 ±0.10 


Peterson et al. (1985) 


0.114 


9 cc;+0.09 


Peterson et al. (1985) 


0.100 


9 t;c;+0.14 


Peterson et al. (1985) 


1.0 - 0.05 


2.735 ±0.06 


Mather et al. (1990) (COBE90) 



Table 1: Recent measurements of the CMBR. 

ii. A hot IGM would Compton scatter photons off of energetic electrons and thereby 
raise the photon energy. 

It should be noted that in both scenarios above, the parameters of the model have been 
adjusted in order to provide a "best fit" for the observations. The first model has great 
difficulty in providing the necessary energy and amount of dust to fit the observed excess. 
The second model is easier to formulate and match to the observations. However this model 
requires that all of the baryonic mass of the Universe be contained in the IGM, eg. Taylor 
and Wright (1989) (hereafter TW) use Qjgm = 0.41, with Hq = 50 km s"'^ Mpc~^. This 
is in direct conflict with the constraint on ^baryon by primordial nucleosynthesis, Steigman 
(1989), Ofo/iQ < 0.02. The claim made by TW that primordial nucleosynthesis constraints 
are "not in strong conflict with our results since extrapolation from current D (deuterium) 
abundance measurements back to primordial abundances is somewhat uncertain since it 
depends sensitively on the amount of stellar processing, and therefore on galactic and stellar 
evolution models." 
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Of course, results from the Cosmic Background Explorer (COBE), Mather et al. (1990), 
hereafter COBE90, seem to indicate that there is no excess in the Wien Tail of the CMBR. 
To quote one astronomer from the AAS meeting of January 1990 (when preliminary COBE 
results were announced), "There is no hot inter-galactic medium." COBE90 is a bit less 
definitive stating only that the hot smooth IGM suggested to explain the cosmic X-ray 
background can be ruled out. The COBE results that support this claim originate from 
nine minutes of observation near the North Galactic Pole by the Far Infrared Absolute 
Spectrometer, FIRAS. The spectral measurements, from 1 cm to 0.5 mm, show less than a 
one percent deviation from an ideal blackbody with equivalent temperature 2.735 it 0.06 K. 
Although the COBE data certainly refutes the findings of the Nagoya-Berkeley Group, 
Matsumoto et al. (1988), the limited wavelength range of observations do not conclusively 
rule out the existence of a hot IGM. 



2. Implications of Primordial Nucleosynthesis for the Physical Conditions in 

Lyman-a Forest Clouds 

The hot big bang model of cosmology, the "standard" model, is now recognized as the 
most accurate description of the Universe's beginning and evolution. Since George Gamow 
and his collaborators put forth the idea of a hot "Big Bang" nucleosynthesis model in the 
late 1940s (Gamow (1946); Alpher et al. (1948); and Alpher and Herman (1950)), many 
authors have extensively studied primordial nucleosynthesis. Appreciable amounts of D and 
He^ are not formed in stellar nucleosynthesis and are destroyed by thermonuclear reactions 
at temperatures above 10^ and 10'' K, respectively. But D and He^ are observable and 
their abundances cannot be explained by non-exotic methods. However conditions are 
present during the first few minutes of a hot expanding Universe to create D and He^. 
Nucleosynthesis began when the temperature of the primitive Universe had cooled to about 
10^ K. This occurred approximately 100 seconds after the Big Bang. Protons fused with 
neutrons to form deuterium nuclei and deuterium nuclei could then fuse to form helium. 
Most of the helium in the Universe was formed at this time, along with deuterium and 
lithium. As the Universe continued to expand and cool, the temperature fell below that 
which is needed to sustain these reactions, and the fusion of D and He^ stopped. 

The number density of the fusion by-products and the rates at which various fusion reactions 
proceed must be known in order to set the limits on when Big Bang nucleosynthesis began 
and ended. This is not a difficult task as these rates are measured in the laboratory. Then 
the evolution of a gas of neutrons and protons, i.e. baryons, is followed in time to calculate 
the results of Big Bang nucleosynthesis. Before nucleosynthesis began, the protons and 
neutrons were in equilibrium and were present in equal numbers. Neutrons could change 
into protons and vice versa via collisions with members of the lepton family, electrons and 
neutrinos and their associated anti-particles. The neutron to proton rate is slightly larger 
than the proton to neutron rate because the neutron is slightly heavier than the proton. 
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However, early in the Big Bang, these reactions proceed at nearly the same rate because of 
the high temperature. 

The ratio of neutrons to protons fell from one to less than a third as the temperature fell to 
10^ K (because of the different reaction rates). At a temperature of 10^ K, this ratio fell 
below a seventh, and the temperature was sufficiently low to allow fusion to begin. The first 
fusion product was deuterium, one proton and one neutron in the nucleus. Interaction of 
deuterium with other protons and neutrons formed tritium (one proton and two neutrons) 
and He^ (two protons and one neutron). These nuclei continued to interact to form He^. 
He^ is a very stable and tightly bound nucleus, consequently almost all of the neutrons were 
bound into He^ nuclei, with traces left in D, He^, Be'' and Li"^. Nucleosynthesis ceased at 
He^ since no stable nuclei is formed from the interaction of He^ and a proton or a neutron. 

The synthesis of these light elements depends only on one adjustable parameter, the density 
of nucleons. A higher nucleon density would "burn" D at a faster rate and hence the D 
abundance would be lowered. Deuterium is only destroyed in the course of galactic chemical 
evolution, therefore an upper limit can be set for the nucleon density from the local observed 
D abundance, Steigman (1989). D is fused by stars into He^, which some survives stellar 
processing. In this manner, the local He^ and D abundances can be used to set a lower 
bound on the nucleon density. These bounds can be refined by using the other light elements 
formed during primordial nucleosynthesis. Steigman (1989) finds that the present density 
of nucleons is bounded by 

0.011 < nbhloo < 0.020 , (2) 

where ilfe is the ratio of baryon mass to the critical mass of the Universe and /i^qq is the 
Hubble constant in units of 100 km Mpc~^. 



3. Properties of Lyman-a Clouds 

The "Lyman-a forest" of absorption lines seen in the spectra of high redshift Quasi-Stellar 
Objects (QSOs) have been the subject of much investigation in recent years. Sargent et al. 
(1980) concluded that this forest of absorption lines resulted from discrete clouds of ionized 
gas along the line-of-sight and are not associated with the QSOs. Therefore these cloud 
systems can be used as windows for examining the Inter-Galactic-Medium as well as other 
cosmologically significant properties. 

Typical physical characteristics for a Lyman-a cloud are: size, ~ 15 kpc; temperature, 
Tg ~ 40,000 K; and neutral hydrogen column density, N^i ~ 10^^ cm~^. The clouds 
are assumed to be photoionized by the cosmic background with a total hydrogen density 
n^* ~ 10~^ — 10~^ cm~^ and in pressure equilibrium with a confining hot phase. Several 
authors have also noted that the number of clouds per unit redshift increase with redshift 
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/? 


7 


Reference 


1.89 ±0.14 


1.7 ± 1.0 


Atwood et al. (1985) 


1.61 ±0.1 




Tytler (1987) 


1.68 ±0.10 




Carswell et al. (1984) 


1.71 




Carswell et al. (1987) 


1.65 ±0.10 




Hunstead et al. (1987) 




2.36 ± 0.36 


Peterson (1983b) 




2.20 ±0.40 


Peterson (1983a) 




2.17 ±0.36 


Murdoch et al. (1986b) 




1.40 ±0.70 


Carswell et al. (1982) 



Table 2: Recent measurements of /3 and 7. 

and decrease with column density described by a distribution function of the form 

P{Nhi, z)dNHidz = BN~1{1 + zY'dNnidz , (3) 

where /?, 7 and B are constants to be evaluated from observations. The expression is valid 
over a ranee of column densities 10 cm^ ~^ < Nhi < W^^ cm-2, Hunstead (1988). Many 
authors have measured /? and 7 as shown in Table [2j The constant B is then evaluated 
from observations and shown in Table [3l This work uses /3 = 1.7±0.2, 7 = 2.3±0.4 and 
B = 1.53 X lO^*' and all calculations will be done at redshift z = 2.5. Therefore Eq. ([3]) can 
be written as 

P{Nhi, z)dNHidz = 1.53 x 10^°iV^}-^(l ± zf-^dNuidz , (4) 
and the total neutral hydrogen density per unit redshift is 

^^Nhi = j NhiP{Nhi, z)dNHi = 5 x lO^^^^s cm'^ , (5) 

where is the uncertainty in the measured neutral hydrogen column density in units of 
5 X 10^*^ cm-2. 

Observations cannot reveal the ionized column density, Nhii, through a cloud, only the 
neutral column density, Nfji- However, the former is needed to estimate the mass of a 
cloud since the cloud is mostly ionized. The ionized and neutral hydrogen number densities 
are related through a simple ionization balance equation 

UeUpaiTe) = T{z)nHi , (6) 

where a and T are the recombination coefficient and photoionization rate, respectively. 
The clouds are much larger than the angular diameter of the QSO. Therefore the average 
chord length through a cloud along the line-of-sight is {D) = D/^/2. This suggests particle 
densities uhi = 3.05 x 1Q~^D^^ cm~^ at z = 2.5, and D^^ is the uncertainty in cloud 
diameter in units of 10^^ kpc. The ionization balance equation is normalized by assuming 
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B X IQio 


Method 


Reference 




1.53 


Total (high z) 


Carsweh (1988), Fig. 


4 


0.63 


Statistical 


Carswell (1988), Fig. 


4 


3.16 


Total (low z) 


Carsweh (1988), Fig. 


4 


0.01 


Statistical 


Carsweh (1988), Fig. 


4 


2.17 


Total 


Hunstead (1988), Fig. 


, 4 


0.08 


Statistical 


Hunstead (1988), Fig. 


, 4 


0.17 


Total 


Atwood et al. (1985) 




0.398 


Statistical 


Tytler (1987) 




10.6 


Statistical 


Carsweh et al. (1984) 




0.855 


Total 


Carswell et al. (1987) 




0.128 


Statistical 


Carswell et al. (1987) 





Table 3: Determination of the constant B. Total method is normalized to neutral hydrogen 
column density and statistical method is normalized to linear least squares regression. 



the clouds have primordial abundances, 90% hydrogen and 10% helium by number; helium 
is fully ionized, Ug = 1.2np; and a typical temperature of Tg = 40, 000 K, with corresponding 
Case B recombination coefficient, aB{Te) = 7.7 x 10"^"* cm^ s~^. The photoionization rate 
is defined as 

f°° T 

'<^' = L ft^'""" • 

Bajtlik et al. (1988) solve for the background intensity from the proximity effect found in 
Lyman-a forest lines. Their measured intensity is log J,^^ = —21.0 it 0.5, where Jj^'^ is the 
Lyman- limit intensity in units of ergs cm~^ s~^ Hz~^. Using this form for the intensity, Eq. 
(HI) is then 

T{z) = 2.99 X 10~^2 J21 s"^ , (8) 

where J21 is the uncertainty in the Lyman-limit diffuse background in units of 47r x 
10~^^ ergs cm~^ s~^ Hz~^. Eq. ([6]) can then be solved for the proton number density 

Hp = 9.93 X 10""^ Jg/^ D'^^"^ cm-3 . (9) 

Assuming the clouds retain their integrity, expand isothermally with the Universe and are 
in pressure equilibrium with a hot phase which is nonrelativistic (at z = 2.5) and expanding 
adiabatically, i.e. P{z) = ntot{z)T{z) ~ (1 + z)^, the electron to neutral hydrogen density 
ratio is 

3.91 X 10^ D~^''^ . (10) 



The electron column density per unit redshift is then 

-Ne = —NhI- 



dz dz uhi 

= 1.96 X 10^1 D~^'^ A5 cm-2 . (11) 
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4. Baryonic Mass Budget of the Universe at z = 2.5 

The mass density contained within the Lyman-a clouds can be determined from Eq. (|lip . 
In an Einstein-de Sitter Universe 

Since the clouds are of primordial abundances, the number density of any one cloud is then 

Pcloud = nipUH + man He , 

= 1.2mpne . (13) 

The mass density of all of the clouds at redshift z = 2.5 is then the mass of one cloud x 
the number density of clouds 



dN^ dz 

p,louds = 1.2m, — - . (14) 



This is evaluated to be 

Pciouds = 9.72 X 10-30 /iioo 4/' D-^''' g cuT^ . (15) 
This result was used to determine the baryonic mass budget of the Lyman-a clouds 

n^iouds = = om2h^,UU'Dli'A, . (16) 

Pcritical 

This result can be used to determine the conditions in the Inter-Galactic Medium at the 
epoch z = 2.5. 



5. Density and Temperature of the IGM at the Epoch z = 2.5 

The baryonic mass budget of the Universe remains an open question. Studies of primordial 
nucleo synthesis, Steigman (1989), show that the fraction of the critical density contained 
in baryons is 

= 0.02 /i^o^ ^ 

where ftioo is the uncertainty in the Hubble parameter, Hq = 100 x hiQQ km Mpc"^. 
It would be surprising if galaxy formation were 100% efficient, and in fact stars and other 
luminous material in galaxies contribute r^gaiaxios ~ 0.007 — 0.017 at the present epoch. Tang 
et al. (1986) and Faber k Gallagher (1979), i.e. 19% to 47% of 0^, for /iioo = 0.75. This 
suggests that the majority of the baryonic mass might occur in another form, perhaps in 
the IGM. Results from the previous section indicate that the baryonic mass contained in 
the Lyman-a clouds is 

^^clouds = 0.6/iioo^^6 • (18) 
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The Lyman-a clouds have been extensively studied. Peterson (1978) and Sargent et al. 
(1980) have shown that these clouds: are not associated with galaxies as evidenced by 
the lack of clustering, have low metal abundance and exist in large numbers. They are 
inter-galactic in nature and of primordial abundances. They are thought to be confined by 
pressure equilibrium with the rest of the IGM, niQM^iGM = f^cioadTdouA- We now assume 
that the remainder of the baryonic mass, as indicated by primordial nucleosynthesis, is 
contained in the rest of the IGM, fi;, = rigaiaxics + ^clouds + ^^igm- Since the IGM is also of 
primordial abundance, the number density of the IGM is 

^^IGM = 0^(1 - 0.35/i2oo - O.eO/iioo) , (19) 

PlGM = Pcritical X ^^IGM j (20) 
PIGM .^-1 N 

""^^^ = 12^ ■ ^^^^ 

, ^ I I tip 

The above equations can be evaluated for niGM 

niGM = 7.42 X 10^^(1 - 0.35/i?oo - 0.60/iioo) ■ (22) 



A lower limit can be set on the temperature of the IGM, at z = 2.5, by choosing the 
minimum mass in the luminous matter, 0,g and assuming pressure equilibrium with the 
Lyman-a clouds 

5.35 X 106 

^^^'^ - 1 - 0.35/.?oo - 0.60/.100 ■ ^ ^ 

With /iioO) ?IGM > 1.5 X 10'' K at z = 2.5. The lack of observed neutral hydrogen, Gunn 
&; Peterson (1965), indicates that if the IGM exists, beyond Lyman-a systems, it is ionized 
and hot as indicated by the above temperature relationship. 



6. Thermal History of the IGM 



We now seek solutions to the temperature-redshift relation for the IGM. The temperature of 
the IGM at any epoch will directly affect the interaction of the cosmic microwave background 
radiation (CMBR) with the hot electrons of the IGM. We have already established the 
temperature of the IGM at the epoch z = 2.5, and this will be used as a boundary condition 
for the temperature-redshift relation as a function of redshift. 

Consider the Universe as a closed system and any expansion or contraction is consequently 
adiabatic. Therefore solutions with constant entropy are desired for the thermal history of 
the IGM. The entropy per particle of a non-relativistic ideal gas is 

| = ^lnr + lny , (24) 

where k is Boltzman's constant, T is the temperature of the gas and V is the volume of 
the gas. Non-relativistic expansion indicates that T ~ and since the volume scales 
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as V ~ (1 + z) ^, the temperature scales as T (1 + -2)^, and the case of coohng a 
non-relativistic gas by adiabatic expansion is solved, Field & Perrenod (1977), hereafter 
FP. 

Unfortunately for this result, the IGM electrons are already nearly relativistic at z = 2.5 
as determined previously, and will become relativistic as the Universe contracts and the 
temperature increases. Therefore the cooling of a relativistic gas by adiabatic expansion 
must also be considered. A dimensionless temperature, Q = kT/rrif.c^ is defined such that 
9 = 1 for T = 5.94 x 10^ K. The general solution for the entropy per particle of any gas is 
given by 

s f^dEde ^ 

+ \nV , (25) 



k J So do e 

where the energy is expressed in units of mgC^. For a purely relativistic gas, E = 36, dE/dO 
= 3, and the entropy per particle is 

f=3.n(^)+l„P' . (26) 

Therefore in an adiabatic expansion, 6 ^ V~^l^ ~ (1 + z). The IGM consists of electrons , 
protons and a-particles. The protons and a-particles will always be non-relativistic for the 
temperatures considered here, but the electrons may become relativistic. The entropy per 
particle for this gaseous mixture is then 



S rdEdO ^ 10 



?ln( — 1 + lny 



(27) 



In the non-relativistic limit, E = 3kT/2, and in the relativistic limit, E = 3kT. The 
dependence of the electron energy upon temperature in the transition region between these 
two limits must be known to solve Eq. (j25p . 



Consider any quantity yl of a system of particles, of mass m, that is dependent upon the 
momentum of the particles. Then the average value of A, (A), Tolman (1917), is 

{A)=4ttvJ^ aexp (^VV'^ + m'^c'^) Aip'^dip , (28) 

where a is a constant and ip is the momentum. The energy, in units of mgC^ is then 

{E) = ^ r aexp ( -^V'02 + m2c2') ^^2 + ^2c2^2^^ ^ (29) 
mc Jq \kT J 

The unknown constant a may be eliminated by normalization 

47rF /" aexp (y^VV'^ + Tr{^cA i?di) = 1 . (30) 

The substitution 



7=\/l + 4^ , (31) 



simplifies Eq. (j29p and the average kinetic energy per particle of the electron gas becomes 



which corresponds to Eqs. (8) and (9) of TW. TW solved the denominator of Eq. ()32p 
by integrating in parts to achieve 6K2{l/9), where K2 is a modified Bessel function of the 
second kind. To complete the integral in the numerator of Eq. (j32p . we note from the 
recursion relations for modified Bessel functions, Abramowitz &: Stegun (1964), 

^K,{z) = ^K2{z)-K3{z) . (33) 
Integration by parts with the above relation, the numerator of Eq. (132p is 



e-'^/^Vl^ - Idj = OKsil/O) - e^K2{l/e) , (34) 
and the average kinetic energy is 



{E) = e 



e \K2{i/e) 



a{9)e . (35) 



The function a{6) is shown in Figure[TJ The derivative of energy with respect to temperature 
(needed for the entropy per particle calculation) is 

dE _ 1 / j^3(i/g) y K^ji/e) 

de \eK2{i/e)) ^9x2(1/9) ' ^ ' 

The volume may be eliminated from Eq. ([25]) by 1/ = Vq{1 + z)-'^. Then Eq. ^ is used to 
integrate Eq. ()25p to define the entropy per particle for both relativistic and non-relativistic 
adiabatic cooling by expansion. An additional source of cooling for the electron gas in the 
IGM is from inverse Compton scattering of CMBR photons off of hot IGM electrons. The 
change in entropy with respect to redshift from inverse Compton scattering is 

ldSc^^J_dQ^ld^ ^^^^ 
k dz kT dz 9 dz 

The Compton power loss per scattering, Rybicki & Lightman (1979), is 



dE 4 2n2 Uph 

at 6 mpC^ 



CTTCT^/?'^ , (38) 



where ax is the Thompson scattering cross section, /3 has its usual relativistic meanings, 7 is 
as defined before, and Uph is the initial photon energy density, in this case the undisturbed 
CMBR, Uph = aT^, where is the photon temperature. Therefore we can write the 
entropy change for one electron due to inverse Compton scattering as 
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0.001 



1000 



Fig. 1. — The function a{9). Note that at low 6, i.e. non-relativistic temperatures, a = 3/2 
and the average kinetic energy is 3kT/2. For large 9, i.e. relativistic temperatures, q = 3 
and the average kinetic energy is 3kT. The transition between the two limits is smooth, 
"nature abhors sharp corners," Professor of Physics Emeritus Stan Christensen, Kent State 
University. 



Setting = 1 implies dt/dz = -Hq^{1 + z)~^/^. The power loss per electron is averaged 
over the electron distribution 



1 dS, 



1 



k dz 



(40) 



where the substitution = r^o(l + -2), as the CMBR photons are always relativistic. The 
integral is evaluated to be 



k dz Horuec^ ^ K2{l/e) 



(41) 



The total entropy per particle of the IGM is then 



S_ 

k 



dEde 15 /e\ 63, 
-TTT^ H 111 111 



do 



do e 



11 



11 



l + z 

1 + ^0 



+ 



1 dSr 



k dz 



-dz 



(42) 



The temperature 9 required to hold Eq. (j42p constant is found as a function of redshift. This 
defines the temperature-redshift relationship for the IGM, Figure [2j Also shown in Figure 
[2] is the case for non-relativistic expansion and the model calculated by TW. The boundary 
conditions in Figure [2] are 9 = 0.00252 at z = 2.5, as found previously (lower limit), and 
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Temperature-Redshift Relation for the IGM 



Taylor Wright (1989) 




I I n < U.--.----I I I I I 

2 4 6 8 10 12 14 16 

redshift (z) 

Fig. 2. — The thermal history of the IGM for Hq = 75 km Mpc~^, for different initial 
temperatures 9. Also shown is the model proposed by TW and the temperature from purely 
non-relativistic adiabatic expansion. 

9 = 0.04 as an intermediate value. The difference between this work and TW is evident 
in Figure [2l The differences result from the boundary conditions. TW chose 9q in order 
to fit the observed x ray background. This work calculates from pressure equilibrium 
between the Lyman-a cluds and the constraints on the fractional mass contained in baryons 
by primordial nucleosynthesis. It is clearly evident that the model proposed by TW would 
have a substantially greater effect upon the CMBR than the model presented here. 



7. Distortion of the CMBR by the Hot IGM 

We have previously shown that the IGM, if it exists, is hot and optically thin. Inverse 
Compton scattering, an efficient cooling mechanism for the IGM, also redistributes energy 
from the hot electrons to the photons of the CMBR. The redistribution of energy will in 
turn change the distribution in frequency of the photon gas, hence distorting the blackbody 
spectrum of the relic radiation. The effect is most noticeable in the Wien tail of the spec- 
trum, and many authors have previously calculated this distortion in efforts to understand 
the reported excess in the Wien tail of the CMBR. Therefore this work is not unique in 
method, but because of the vastly different thermal history of the IGM presented in the 
previous section as compared to other author's efforts (most notable TW and FP), our 
calculated distorted CMBR is correspondingly different, albeit of the same nature as TW 
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and FP. 

The reader is directed to TW for a superb description of the methodofogy for calculating the 
distortion of the CMBR, it will not be reproduced here. Only the results will be reported. 

The results are shown in Figure [3l Three different models are shown in Figure [3] correspond- 
ing to different epochs when radiation from QSOs first ionized the IGM. These models are 
Model 1 - Zon = 6, Model 2 - Zon = 15, and Model 3 - Zon = 20. The models start 
by adjusting the initial temperature of the CMBR so that the output temperature of the 
distorted CMBR is T^q = 2.735 K in the range 1 cm > u-^ > 0.05 cm, to correlate to 
COBE90. (This is done only to show that the results of COBE90 do not disallow a hot 
IGM.) Previous authors have adjusted the density and temperature of the IGM in order 
to model the cosmic x ray background and the highly distorted CMBR by Matsumoto et 
al. (1988). The initial scope of this work was to adjust the parameter Zan in order to 
present a model of the IGM that also matched the results of previous measurements of the 
CMBR, discounting COBE90. The preliminary results of COBE90 have changed the path 
of this work. Whereas it has been boldly stated that COBE90 has shown there is not a 
hot (ionized) IGM, the results shown here show that COBE90 data is inconclusive as to 
the existence of the hot IGM. The COBE90 results do put limits on the parameter space 
of temperature and density of the IGM as evidence by Figure El A hot, thin IGM may still 
exist, but certainly not to the extent as to be the source of the cosmic x ray background as 
has been previously postulated, TW and FP. 



8. Summary 

Previous efforts. Fisher (preprint), calculated an effective recombination coefficient for opti- 
cally thin hydrogen clouds. An evolutionary model of the diffuse UV background was then 
calculated and its effects on the stability of metal-free Lyman-a cloud systems, because 
of radiation pressure, was modeled. From this model, limits were derived on the epoch of 
photoionization of the Lyman-a clouds. Primordial nucleosynthesis and the observed mass 
of luminous matter in the Universe were used to determine the temperature and density 
of the IGM at the epoch z = 2.5. This was used as a boundary condition to determine 
the thermal history of the IGM. Finally, the distortion of the CMBR because of this hot 
IGM was calculated and compared to previously published data. This mode was found 
to disagree with the distortion as measured by the Berkeley-Nagoya group, but was found 
to be consistent with the preliminary results of the Cosmic Background Explorer to high 
accuracy. This model of the IGM is not capable of producing the cosmic x ray background 
as measured. 

The most interesting result of this work is its implication for the distortion of the CMBR. 
Much effort had been made to explain the distortion in the Wien tail of the CMBR, as 
measured prior to the launch of the COBE satellite. The original intent of this study was 
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Fig. 3. — Three different models for the distortion of the CMBR by inverse Compton 
scattering off of hot IGM electrons is shown above. All three models have 6{z = 2.5) = 
0.0025. The models differ in the epoch at which the radiation from QSOs first ionized the 
IGM, these are Model 1 - Zon = 6, Model 2 - Zon = 15, and Model 3 - Zon = 20. Also 
shown are previous temperature measurements of the CMBR by previous authors including 
COBE90. It is noted that all three models from this work are in agreement with the 
COBE90 measurements (including the uncertainty in the measurements). 



to proceed along the same path as previous authors, most notably TW, and determine 
if the then observed distortion could be attributed to inverse Compton scattering of the 
CMBR photons by hot IGM electrons. The hoped for end result was to compare the 
required density and temperature of the IGM to the observed conditions in the Lyman-a 
forest cloud systems. Only after the model was complete, and the lack of a gross distortion 
similar to that published by the Berkeley-Nagoya group failed to appear, were the results 
of COBE made known. Whereas it has been stated that COBE disproves a hot IGM, this 
study indicates that the results, published to date, only dispel the previous results of the 
Berkeley-Nagoya group, and set a boundary, in temperature-density space, allowed for the 
IGM. A hot IGM may still exist and only after COBE extends the spectra to even shorter 
wavelengths, ~ 0.01 cm, will the limits on the existence of the hot IGM proposed here be 
discounted or substantiated. 

Several observations are needed to support the theoretical findings presented here, most 
notably the COBE spectra noted above. Other observations are needed to better define the 
limits on many of the physical characteristics of the Lyman-a cloud systems, including size. 
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temperature, and distribution in neutral hydrogen column density. All of these values were 
used to determine the temperature and density of the IGM. However, it is noted that drastic 
changes in the values of these parameters is not expected, and therefore the result of a hot 
IGM will not change. The results of extending the COBE spectra to shorter wavelengths 
is most anxiously awaited. 

The author gratefully acknowledges J. A. Baldwin, E.R. Capriotti, G.J. Ferland, J.V. Vil- 
lumsen and R.J. Weymann for helpful dscussions and comments on the manuscript. This 
work was supported in part by National Science Foundation Grant AST 000-00. 
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